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We investigate the virtual photon structure function in the supersymmetric QCD (SQCD), where we have squarks and gluinos in 
addition to the quarks and gluons. Taking into account the heavy particle mass effects to the leading order in QCD and SQCD we 
evaluate the photon structure function and numerically study its behavior for the QCD and SQCD cases. 
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Since the experiments at the Large Hadron Collider (LHC) 
(jl started there has been much anticipation for the signals of 
the Higgs boson as well as for an evidence of the new physics 
beyond Standard Model such as supersymmetry (SUSY). Once 
these signals are observed more precise measurement needs to 
be carried out at the future e + e~ collider, so called International 
Linear Collider (ILC) J2l. In such a case, it is important to know 
the theoretical predictions at high energies based on QCD. 

It is well known that, in <? + e~ collision experiments, the cross 
section for the two-photon processes e + e~ — > e + e~ + hadrons 
dominates at high energies over the one-photon annihilation 
process 03J. We consider here the two-photon processes in the 
double-tag events where both of the outgoing e + and e~ are de- 
tected. Especially, the case in which one of the virtual pho- 
ton is far off-shell (large Q 2 = -q 2 ), while the other is close 
to the mass-shell (small P 2 - -p 2 ), with A 2 <k P 2 «: Q 2 
(A: QCD scale parameter), can be viewed as a deep-inelastic 
scattering where the target is a virtual photon and we can cal- 
culate the photon structure functions in perturbation theories 
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Some time ago the effects of supers ymm etry on two-photon 
process were studied in the literature lil3l[l4 , 15 . Hill . In this 
paper based on the framework of treating heavy parton distri- 
butions IB [l8ll we reexamine the effects of the squarks and 
gluinos appearing in SUSY QCD (SQCD) on the photon struc- 
ture functions to the leading order in SQCD which can be mea- 
sured in the two-photon processes of e + e~ collision illustrated 
in FigfTJ 



1. Evolution equations for the SUSY QCD 

We consider the DGLAP type evolution equations for the 
parton distribution functions inside the virtual photon with the 
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Figure 1: e + e~ two-photon processes in supersymmetric QCD. The solid 
(dashed) line denotes the quark (squark), while the spiral (spiral-straight) line 
implies the gluon (gluino). 



mass squared, P 2 , in SQCD where we have squarks and gluinos 
in addition to the ordinary quarks and gluons. Evolution equa- 
tion to the leading order (LO) in SQCD reads as in QCD 

dqy(t) 



q y (t) ® P 



(0) 



r,(0) 



(1) 



dt ' a s (i) 

where and k®' are 1-loop parton-parton and photon-parton 
splitting functions, respectively (see Appendix). The symbol ® 
denotes the convolution between the splitting function and the 
parton distribution function. The variable t is defined in terms 
of the running coupling a s as [20]: 

a s {P 2 ) 
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with the parton distributions probed by the virtual photon with 
mass squared Q 2 as 



q?(t) = (G,A,qi 



(3) 



■q nj , s\ 

where rif is the number of active flavors. In eq.(f2|, fio = 9 - rif 
for SQCD. We denote the distribution function of the i-th flavor 
quark, squark by qi(x,Q 2 ,P 2 ), Si(x,Q 2 ,P 2 ), (i = l,--- ,n f ), 
and the gluon, gluino by G(x, Q 2 , P 2 ), A(x, Q 2 , P 2 ), respec- 
tively. The 1-loop splitting functions were obtained in l2ll l22ll . 
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We first consider the case where all the particles are massless. 
Although this is an unrealistic case, it is instructive to consider 
the massless case for the later treatment of the realistic case 
with the heavy mass effects. 

For the massless partons the evolution starts at Q 2 = P 2 and 
hence we have the initial condition q y (t — 0) = 111 ill . 

The 1-loop splitting function is given by (see Appendix A) 
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where Pab is a splitting function of B parton to A parton with 
A,B = G, A,q and s. While the splitting functions of the photon 
into the partons G, A, q and s, are denoted as (see Appendix B) 



k (0) = (k G ,k A ,k q ,k s ) 



(5) 



We introduce the flavor-nonsinglet (NS) combinations of the 
quark and squark distribution functions as 



n f 



q NS (x, q\ p 1 ) = Ys^i - <e 2 »<?,■(*, e 2 , p 2 ) . (6) 

!=1 

"/ 

s NS (x, Q\P 2 ) = J](e 2 - (e 2 )) Si (x, Q 2 ,P 2 ) , (7) 



!=1 



where e-, is the z'-th flavor charge and (e 2 ) = Yii e2 l n f is tne 
average charge squared. We also define the flavor-singlet (S) 
combinations for quarks and squarks 



X(x,Q 2 ,P 2 ) = Y j q i (x,Q 2 ,P 2 ) 

s(x,Q 2 ,p 2 ) = Y, s t x > e 2 ^ 2 ) 



(8) 
(9) 



We now rearrange the parton components of q y (t) using the 
above flavor non-singlet and singlet combinations as: 



q y (t) = (G,A,1,S,q NS ,sm) ■ 
Then we have the following splitting function 



(10) 



p(0) _ 



' Pec 


Pag 


PqG 


PsC 


Pga 


Paa 


PqA 


PsA 


Pcq 


?Aq 


P 


P 

1 sq 


P Gs 


Pa* 


p 

1 qs 


P 

1 ss 











m 

qs 



l sq 
Pss ) 



(11) 



Thus for the flavor-nonsinglet parton distributions 



9ns ~ (9ns . s ns ) > 



satisfy the following evolution equation: 



d 9l 



dt 



NS = y 
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where the splitting functions are 
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For the flavor-singlet parton distribution 

q y s = (G,A,Z,S) , 
we have 

d ~f = 9>Pf + 4t 

dt s s a s (t) 

where 
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= (Kf\ K?>), Kf> = ^ k qt , K ( p = J] k sr (19) 
i=i !=i 
Now we should notice that there exist the following super- 
symmetric relations for the splitting functions 112 ill : 

(20) 
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Hence if we introduce the following combinations 
<p = Y. + S, V = G + A, 



(24) 



then we obtain the following compact form for the mixing of 
the flavor-singlet part: 
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where we denote K\ — K q + K) and the flavor-nonsinglet 
part <Pns = qNS + ^ns becomes 
d(pNS n . a 



dt 
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where we introduced K 
flavor singlet and non-singlet parton distribution functions we 
can express the virtual photon structure function F y 2 as 

F y 2 (x, Q 2 , P 2 ) = xJ] e 2 ( qi (x, Q 2 , P 2 ) + Si (x, Q 2 , P 2 )) 

i 

= x Yj,e 2 - <e 2 ))(qi(x, Q 2 , P 2 ) + Si (x, Q 2 , P 2 )) 

i 

+x{e 2 ) J^iqtix, Q 2 , P 2 ) + Si (x, Q 2 , P 2 )) 



x(q N s (x, Q 2 ,P 2 ) + sns (x, Q 2 , P 2 )) 
+x(e 2 ){Z{x, Q 2 , P 2 ) + S (x, Q 2 , P 2 )) 
xcpNsix, Q 2 , P 2 ) + x(e 2 )4>(x, Q 2 , P 2 ) 



(28) 



In Fig. 2 we have plotted the virtual photon structure function 
F y (x, Q 2 , P 2 ) in the SQCD as well as in the ordinary QCD for 
Q 2 = (1000) 2 GeV 2 and P 2 = (10) 2 GeV 2 . We have also shown 
the quark as well as the squark components of the virtual photon 
structure function FZ in the case of the SQCD. In contrast to 
the QCD, the momentum fraction carried by the quarks in the 
SQCD case decreases due to the emission of the squarks and the 
gluinos. Hence the jt-distribution of the quarks for the SQCD 
increases at small-* and decreases at large x, i.e. it becomes 
more flat compared to the QCD case as seen from the Fig. 2. 
Adding the two components together we get the FZ structure 
function for the SQCD which shows a behavior quite different 
from that of the QCD. 

Q 2 =(1000) 2 GeV 2 , P 2 =(10) 2 GeV 2 



The evolution equations and the splitting functions read 
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Figure 2: The virtual photon structure function F^x, Q 2 ,P 2 ) divided by the 
QED coupling constant a for massless QCD (solid line) and SQCD (dashed 
line) with n f = 6, Q 2 = (1000) 2 GeV 2 and P 2 = (10) 2 GeV 2 . Also shown are 
the quark (dash-dotted line) and the squark (double-dotted line) components. 

2. Heavy parton mass effects 

Many authors have studied heavy quark mass effects in the 
nucleon |23] and the photon structure functions |24, 25l EfSll . 
Now we consider the heavy parton mass effects, and we decom- 
pose the parton distributions in the case where we have «/ - 1 
light quarks and one heavy quark flavor which we take the n y-th 
quark and all the squarks have the same heavy mass, while the 
gluino has another heavy mass |17lll8ll : 



q y (t) = (G,A,qw- ,q„ f -i, s u ■ ■ ■ ,s nf -uqH,s H ) ■ (29) 

We denote the z'-th light flavor quark, squark by q,(x, Q 2 , P 2 ), 
Si{x,Q 2 ,P 2 ), (i = 1, ••• ,«/ - 1), one heavy quark and its 
superpartner (squark) by qn, sh and the gluon, gluino by 
G{x, Q 2 , P 2 ), A(x, Q 2 , P 2 ), respectively. 

We now define light flavor-nonsinglet (LNS) and singlet (LS) 
combination of the quark and the squark as follows: 

«f— 1 Hf— 1 

qLNS = ~~ qi, "LNS = 2^ ( e ? ~ ^ 2 )i) S ! ' 
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i=l i=l n f 

Then we rearrange the parton distributions as 
q y (t) = (G, A, q LS , s LS , q H , s H , qLNS , $lns ) 
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While the photon-parton splitting functions are 

= (kc, kx, k qLS , k SLS , k qH ,k SH , k C[LNS , k StNS ) . 



(32) 



(33) 



Now we take into account the heavy mass effects by setting the 
initial conditions for the heavy parton distribution functions as 
discussed in [18, 271 12811 . 



We note here that the structure function F^ can be written 
as a convolution of the parton distribution q 7 (x, Q 2 ,P 2 ) and the 
Wilson coefficient function C(x, Q 2 ): 



F\(x, Q 2 ,P 2 )/x = q y ®C . 
The moments of the parton distributions are defined as 



q r (n,t)= f dxx"- l q y {x,Q 2 ,P 2 ) , 
Jo 



(34) 



(35) 



where we put the initial conditions: 

q y (n, t = 0) = (0, l(n), 0, fo(n), q H (n), s H (n), 0, s ms (n)) , (36) 
and require that the following boundary conditions are satisfied: 
A(n, Q 2 = m\) = 0, s LS (n, Q 2 = m 2 q ) = 0, q H (n, Q 2 = m 2 H ) = 0, 

(37) 



Sh{ti, Q 2 = m 2 sq ) = 0, s LNS (n, Q 2 = m 2 sq ) = 0, 



where m^, m sq and m# are the mass of the gluino, squarks and 
the heavy (here we take top) quark, respectively. Note that here 
we take all the squarks have the same mass m sq . 

By solving the evolution equation taking into account the 
above boundary condition we get for the moment of q y : 



q y (n,t) = 
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where the P" is the projection operator onto the eigenstate A, of 
the anomalous dimension matrices y n : 



(39) 



where the anomalous dimension matrices % is related to the 
splitting function P(x) as 



% = -2 



dxx"- l P(x), 



(40) 



and d" = A" /2fio. is the anomalous dimension correspond- 
ing to the photon-parton splitting function: 

Kf> =2 f dx^ 1 ^*). 
Jo 

The initial value q y (n, 0) is determined so that we have 

2 a s (P 2 ) 
q)(t = t mj ) = 0, t m = — In - , 
P0 a s (mj) 
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for j = A,sls, 1h, $h and Slns ■ By solving the above 
coupled equations we get the initial condition: q y (n,Q) = 
(0, A{n), 0, s LS (n), q H (n), s H (n), 0, s LNS («)). 

Now we write down the moments of the structure function 
in terms of the parton distribution functions and the coefficient 
functions, which are 0(a®) at LO. We take 



Cf ( 1 , 0) r = (0, 0, (e\, (e\,e 2 H ,e 2 H ,l,l). 



(44) 



Then the n-th moment of the structure function F 2 to the lead- 
ing order in SQCD is given by 

Ml = f dxx n - x F y 2 lx = q y (n) ■ C Y „(l,0) 
Jo 

= <<? 2 >l<?ls + {e 2 )iSLS + e 2 H q 2 H + e 2 H s H + q LNS + s LNS ■ (45) 



3. Numerical analysis 

We have solved the equations (l43l for q y (n, 0) numerically, 
and plug them into the master formula ( l38l l for the parton distri- 
bution functions and then evaluate the moments of the structure 

By inverting the Mellin 



function FZ based on the formula 
moment we get the F^ as a function of Bjorken x. 

In Fig. 3, we have plotted our numerical results for the Fl/a. 
The 2dot-dashed and dashed curves correspond to the F y 2 /a 
for the massless QCD and SQCD, respectively, where all the 
quarks and squarks are taken to be massless. Of course this is 
the unrealistic case we discussed in the previous section. For 
the more realistic case, we take rif = 6 and treat the u, d, s, 
c and b to be massless and take the top quark t massive. We 
assume that all the squarks possess the same heavy mass and 
the gluino has another heavy mass. In these analyses, we have 
taken Q 2 = (1000) 2 GeV 2 and P 2 = (10) 2 GeV 2 . For the mass 
values we took the top mass m top - 175 GeV, the common 
squark mass, m s = 300 GeV and the gluino mass nix = 700 
GeV. 

The double-dotted curve shows Fl/a for the QCD with the 
mass of the top quark as well as the threshold effects taken into 
account. The dash-dotted curve shows the quark component 
for the massive SQCD case with massive top quark, while the 
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Figure 3: Fl(x, Q ,P )/a with SUSY particles as well as top threshold in- 
cluded. The dashed (2dot-dashed) curve corresponds to the massless SQCD 
(QCD) case. The double-dotted curve shows the massive QCD case. The dash- 
dotted (dotted) curve corresponds to the quark (squark) component of the mas- 
sive QCD. The solid curve means the F^/a for the massive SQCD. The kink at 
a— 0.89 (0.74) corresponds to the top (squark) threshold. 



dotted curve means the squark component for the same case. 
The sum of these leads to the solid curve which corresponds 
to F y /a for the massive SQCD with massive top and threshold 
effects included. Here, we adopt the prescription for taking into 
account the threshold effects by rescaling the argument of the 
distribution function f(x) as H29fl : 



f(x) -> f(x/x m . ix ) , 



1 



1 + S + ^ 



(46) 



where x max is the maximal value for the Bjorken variable. After 
this substitution the range of x becomes < x < x max . At small 
x, there is no significant difference between massless and mas- 
sive QCD, while there exists a large difference between mass- 
less and massive SQCD. At large x, the significant mass-effects 
exist both for non-SUSY and SUSY QCD. The SQCD case is 
seen to be much suppressed at large x compared to the QCD. 
The squark contribution to the total structure function in mas- 
sive SQCD appears as a broad bump for x < x max . Here of 
course we could set the squark mass larger than 300 GeV, e.g. 
around 1 TeV, as recently reported by the ATLAS/CMS group 
at LHC, for higher values of Q 2 . 

4. Conclusion 

In this paper we have studied the virtual photon structure 
function in the framework of the parton evolution equations for 
the supersymmetric QCD, where we have PDFs for the squarks 
and gluinos in addition to those for the quarks and gluons. 

We considered the heavy parton mass effects for the top 
quark, squarks and gluinos by imposing the boundary condi- 
tions for their PDFs in the framework treating heavy particle 
distribution functions Hill . The PDF for the heavy particle with 
mass squared, m 2 are required to vanish at Q 2 — m 2 . This can be 



4 



translated into the initial condition for the heavy parton PDFs, 
q y (t = 0). Due to the initial condition the solution to the evo- 
lution equation is altered as given by (l38l l. This change leads 
to the heavy mass effects for the PDFs. As we have shown 
in Fig. 3, there is no significant difference in the small-x re- 
gion between QCD and SQCD, while at large x, it turns out 
that there exists a sizable difference between the massive QCD 
and SQCD. When compared to the squark contribution to Fl in 

I — ~~L 

the parton model calculation [30], the squark component in the 
SQCD is suppressed at large x due to the radiative correction. 
We expect that the future linear collider would enable such an 
analysis to be carried out on photon structure functions. 

Appendix A. Anomalous Dimensions for SUSY QCD 

Note that the our convention for the anomalous dimension is 
related to the above splitting function as 

(A.l) 
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The 1-loop anomalous dimensions for SUSY QCD are given by 
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where Cf = 4/3 and = 3 for SQCD. Hence we have the 
following anomalous dimensions for the supersymmetric case: 
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where we have the following replacement: n jy n ^ v — > y^ y . In 
the case of non-supersymmetric QCD we have the following 
anomalous dimensions: 
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Appendix B. Photon-parton mixing anomalous dimen- 
sions 

The photon-parton splitting function can be connected to the 
photon-parton mixing anomalous dimensions given by 
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